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ABSTRACT 

It is well known that non-ideal magnetohydrodynamic effects are important in the dynamics 
of molecular clouds: both ambipolar diffusion and possibly the Hall effect have been identified as 
significant. We present the results of a suite of simulations with a resolution of 512"^ of turbulent 
decay in molecular clouds incorporating a simplified form of both ambipolar diffusion and the 
Hall effect simultaneously. The initial velocity field in the turbulence is varied from being super- 
Alfvenic and hypersonic, through to trans- Alfvenic but still supersonic. 

We find that ambipolar diffusion increases the rate of decay of the turbulence increasing the 



decay from t ° to t 



The Hall effect has virtually no impact in this regard. The power 



spectra of density, velocity and the magnetic field are all affected by the non- ideal terms, being 
steepened significantly when compared with ideal MHD turbulence with exponents. The density 
power spectra components change from ~ 1.4 to ~ 2.1 for the ideal and non-ideal simulations 
respectively, and power spectra of the other variables all show similar modifications when non- 
ideal effects are considered. Again, the dominant source of these changes is ambipolar diffusion 
rather than the Hall effect. There is also a decoupling between the velocity field and the magnetic 
field at short length scales. The Hall effect leads to enhanced magnetic reconnection, and hence 
less power, at short length scales. The dependence of the velocity dispersion on the characteristic 
length scale is studied and found not to be power-law in nature. 
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The role of turbulence in molecular cloud evo- 
lution has been a subject of much study in the 
lit erature (see, for example, the excellent reviews 
of Mac Low fc KlessenI [20o3 lElmegreen fc Scalol 



20041 ). Observations of the propert ies of gas and 



dust motions in molecular clouds (jLarsonl 119811 ) 
suggest that, indeed, turbulence is present. It is 
clear that turbulent motion could influence the 
star formation rate and efficiency as well as the 
initial mass function ( !Elmegreeiilll993 : Klein et al 



[2003). Given all this there are several interesting 
questions which arise: 
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1. What is the source of molecular cloud tur- 
bulence? 

2. How fast does it decay? 

3. How does it affect star formation? 

4. How does it affect the evolution of molecular 
clouds? 

The first two questions are clearly inter-related: if 
the turbulence decays very quickly then we need a 
lot of energy from its source in order to maintain 
it. Indeed, to address the final two questions we 
must first gain insight into the first two. 

In order to study turbulence in molecular 
clouds we must resort to numerical simula- 
tions. Quite simply, there are no satisfactory 
analytic techniques for addressing compressible 
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maKnetohydrodynamic (MHD) turbulence 
Elmegreen fc Scaloll2004l ). 



(,e-g. 



Many authors have performed sophisticated 
numerical simulations in order to investigate 
both the qualit ative nature of MHD t urbulence 
and i ts decay (iMac Low et al. 1998t Mac Low 



199S; Ostriker et al. 12001: Vestuto et al 



2003; 



2007 



Gustaffson et alJ 1 2006- 'Glover fc Mac Low 
Lemaster fc Stona, 2008, ,2009,) . Most of this work 
has been carried out for the case where ideal 
MHD is valid - i.e. on relatively large length 
scales. When smaller length scales are consid- 
ered (e.g. lengths of significantly less than a par- 
sec) ambipolar diffusion becomes non-negligibl e 
in molecular clouds (.Oishi fc Mac Low! 2006|) . 
Some authors dLi et al.l2008HKudoh fc Basul20o'a 



Oishi fc Mac Lo J200g ) have studied driven MHD 
turbulence in the presence of ambipolar diffusion. 
All these authors find that ambipolar diffusion 
produces significant differences in the properties 
of the turbulence. 

ft has been suggested (IWardld liool that al- 
though the Hall resistivity is generally at least an 
order of magnitude lower than the ambipolar re- 
sistivity in molecular clouds, its effect should not 
be ignored due to the qualitative change it in- 
duces in the behavior of the magnetic field. Re- 
searchers working on reconnection and the so- 
lar wind have studied the Hall effect in the con- 
text of turbulence and found that, although the 
overall decay rate appears not to be affected, 
the usual coincidence of the magnetic and veloc- 
ity fiel ds seen in MHD does not occur at small 



scales ( Servidio et al 



20071 : iMininni et aD l2006t 



Matthaeus et al. Almost no work has been 



done on comparing the influences of this effect 
coupled with that of ambipolar diffusion on tur- 
bulence (with the exceptio n of low resolution sim- 
ulations bv Down es fc O'SullivanI l2008l) . In par- 
ticular, to our knowledge the work presented here 
represents the first systematic study of molecular 
cloud turbulence incorporating both the Hall ef- 
fect and ambipolar diffusion simultaneously. 

The main aim of this work is to examine in 
detail the differences between the decay of ideal 
MHD turbulence and that of more realistic non- 
ideal MHD turbulence with a full tensor resistiv- 
ity incorporating the effects of ambipolar diffusion, 
the Hall effect and Ohmic resistivity. This work 
is new in two respects: no previous work has fo- 



cussed on decaying (i.e. un-driven) turbulence in 
the presence of non-ideal terms and, in addition, 
no previous work has addressed the issue of turbu- 
lence in the presence of both ambipolar diffusion 
and the Hall effect simultaneously. This is the first 
of a short series of papers describing a comprehen- 
sive study of multifluid MHD turbulence in a pa- 
rameter regime appropriate to molecular clouds. 

In section[2]we outline the numerical techniques 
used in this work, as well as the initial conditions 
and general set-up for the simulations while in sec- 
tion [3] we describe the methods used to analyze 
the simulation data. In section [4] we present and 
discuss the results of our simulations of turbulent 
decay. Finally, section [5] contains a summary of 
our results. 

2. Numerical method 



In t his w ork we use the code HYDRA ( O'Sullivan fc Downed 

20061 . 2007 ) to integrate the equations of non- 
ideal MHD (see section [2T|) . We assume that the 
molecular cloud material we are simulating can be 
treated as isothermal and that initially the density 
and magnetic field are uniform. For this work we 
assume spatially and temporally constant Ohmic, 
ambipolar and Hall resistivities (see section [^?T|) . 

2.1. Equations and algorithm 

We briefly outline the equations and assump- 
tions in our model h e re but refer the reader to 
O'Sullivan fc Downed (|2006l l2007l) for a compre- 
hensive description of both the full abilities of the 
HYDRA code and the assumptions underlying the 
equations used. 

We assume that the cloud material can be 
treated as weakly ionized. This is clearly valid 
for molecular clouds and allow s us to ignore the 
inert i a of the cha rged species ( Ciolek fc Roberg3 
2002t lFallell2003f) . The equations solved in this 
work are then 
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where p, q, a, \, B and J are the neutral mass 
density, neutral velocity, sound speed, identity ma- 
trix, magnetic field and current density respec- 
tively. The electric field in the frame of the 
fluid, E' , is calculated from the generalized Ohm^ 
law for weakly ionized fluids (e.g. Falld 20031 : 
[crSullivan fc Downesll2006h and is given by 



E' = Ec 



Eu + Ea, 



where 



£^0 = (J -00)00, 
Ea = J X Oh, 

Ea = -(J X Oa) X OA, 



(6) 



(7) 
(8) 
(9) 



using the deflnitions oq — foB, an — fuB, a a — 
JaB, where /o = ^/ro/B, /h = m/B, /a = 
y/vA/B. ro, rn and ta are the Ohmic, Hall and 
ambipolar resistivities respectively. 

In this work these resistivities are kept constant 
in both space and time. We note that, phys- 
ically, they depend on both the magnetic field 
and the densi ty of the various charged species in 
the fluid (e.g. ICiolek fc Robergell2002l : lFallell2003t 
lO'Sullivan fc Downed I2OO6I l2007t ) and hence, in 
reality, do vary in both space and time. Treating 
the resistivities in such a simple way allows us to 
gain a deeper understanding of their influence on 
turbulence without having to consider the compli- 
cating effects of dynamically varying resistivities 
at the same time. As such, and as a first step 
away from the approximation of ideal MHD and a 
single form of magnetic diffusion we believe this is 
an interesting study. Having gained some insight 
into this simplified model a follow-up paper will 
address multifluid MHD turbulence under the in- 
fluence of self-consistently calculated resistivities. 

As noted bylMlidlool) and lO'Sullivan fc Downed 
( 20061 ). the main difficulty with standard numeri- 
cal techniques for integrating equation [3] lies with 
the Hall term. As this term becomes dominant the 
stable time-step go es to z ero. However, O'Sullivan 
fc Downes (|2006l: 12007? ) presented a novel, ex- 
plicit numerical method for integrating this term 
such that the limit on the stable time-step is not 
overly restrictive. We use this "Hall Diffusion 
Scheme" in this work. Of course, all explicitly 
differenced diffusion terms give rise to a stable 
time-step which is proportional to Ax"^ , where Ace 
is the resolution of the simulation. To ameliorate 



this we use standard sub-cycling of the Hall terms 
and super time-stepping to acc elerate the am- 
bipolar d iffusion term s (see Alexiades et al. 1996t 



O'Sullivan fc Downed[2006Ll2007D . 

Equations [1] - [3] are solved using a standard 
shock-capturing, second order, finite volume, con- 
servative s cheme. Equat i onlHi s enforced using the 



method of iDedner et all (|2002l ). The effects of the 



diffusive terms are then incorporated in an opera- 
tor split fashion. 

2.2. Initial conditions 

We examine the decay of supersonic MHD tur- 
bulence in conditions suitable for dense regions of 
molecular clouds. While the simulations presented 
here are, of course, scale-free we present the initial 
conditions used in standard astrophysical units for 
ease of reading. 

The simulations are carried out in a cube of side 
L — 0.2 pc with periodic boundary conditions be- 
ing enforced on all faces. The sound speed is set 
to 0.55kms^^and the initial density is chosen to 
be uniform with a value of lO^cm"^. The mag- 
netic field is also initially uniform in the (1,1,1) 
direction with a magnitude of 1 mG. For these con- 
ditions, suitable conductivities are co = 1 x 10^° 
s~ \ (Th = 10~^s~^ an d (ta = 10~^ s~^ (see figure 
1. IWardle fc Ng||l999l ). We choose these particu- 
lar physical conditions with a view to maximizing 
the influence of the Hall effect in our simulations 
( Wardle fc Ngl 1999f ) . In this way we can use our 
simulations to flnd whether the Hall effect is ever 
likely to be important in molecular cloud turbu- 
lence. 

The initial velocity field is used to instigate the 
turbulence in these simulations. Each component 
of the velocity field is defined to be the sum of 
waves with 16 wave- vectors, each with random am- 
plitude and phase - i.e. 

4 

1.171,71—1 

(10) 

where i = 0,1,2 defines the component (x, y or 
z respectively) of the velocity, A'™" and (/)'™" are 
the random amplitudes and phases and 



fe^-^(l-<5.) 



(11) 
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where Su is the usual Dirac delta function. The 
inclusion of the term in parenthesis in the defini- 
tion of fc' restricts the initial velocity field to be 
solenoidal. Note that we perform all the analysis 
of these simulations at t > 0.2 tc (1 flow crossing 
time) at which time the effects of the precise initial 
conditions used should be negligible. 

The nomenclature for the simulations is xx-ab-c 
where xx denotes the type of physics (e.g. a stan- 
dard molecular cloud run is "mc" , ideal MHD is 
"mhd" etc) , ab is the initial rms Mach number and 
c is the resolution used. The initial root-mean- 
square (rms) of the field is chosen to be either 
Mach 2.5, 5 or 10 depending on the simulation 
in question. These correspond to Alfvenic Mach 
numbers of approximately 0.96, 1.9 and 3.85 re- 
spectively. In addition to the non-ideal MHD sim- 
ulations described we also run 4 further simula- 
tions. The first is an ideal MHD simulation (mhd- 
5-512) which we use for comparison purposes, an- 
other is a pure hydrodynamic simulation (hd-256- 
0.5), and the other two (ambi-5-512 and hall-5- 
512) only incorporate one of ambipolar diffusion 
or the Hall effect, respectively. We use these lat- 
ter simulations to separate out the effects of each 
of these diffusions to better understand the physics 
occurring. Table [T] contains definitions of the var- 
ious simulations used in this work. 

3. Analysis 

In this section we discuss the method of anal- 
ysis of the output of the simulations described in 
section 12.21 The main aim of this paper is to in- 
vestigate the decay rate of supersonic turbulence 
in molecular clouds. Hence, the main analysis car- 
ried out on the simulation results is the calculation 
of the volume-averaged kinetic, magnetic and to- 
tal energy in the simulation as a function of time. 
These quantities are defined as 
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(12b) 
(12c) 



where the angle brackets denote averaging over 
the computational domain. Note that Cb is there- 
fore the difference in the magnetic energy at 
the current time and the magnetic energy at 
t = (see, for example, Vestuto et al. 20031 



Lemaster fc Stond l2009l) . given the assumption 
that no external electromotive force is applied 
and that our boundary conditions are periodic. 
We also calculate the mass- weighted average Mach 
number, defined by 



-I) 



1/2 



(13) 



where a is the sound speed and the mass- weighted 
velocity dispersions, a a, are defined by 



{P) 



1/2 



(14) 



where a is either x, y or z and the angle brackets 
den ote av eragin g over the co mputational domain 



Lemaster fc Stonell2009l) 



[see 

In section [4.31 we present the power spectra for 
the velocity, density and magnetic field for each of 
the 512^ simulations. These spectra are calculated 
by taking the power spectrum in the x, y and z 
directions and then integrating the power for all 
k satisfying k < \k\ < k + dk for each k with 
dk = 1. This gives us some insight into the scale 
of structures being formed by the turbulence for 
the various initial conditions and range of physics 
examined. 

Finally, in section 14.41 we calculate the velocity 
dispersion as a function of length scale, I. For 
these purposes we define the velocity dispersion 
to be 



'^(0 = {< f^lil) >do 



+ < alii) >d O 



+ < alii) >d O 
(15) 



where 



«(0 = {('Z^),-(9a)f}' (16) 



where < • >; indicates an average taken over a 
cube of side I in the simulation domain and < 
• >domain indicates averaging of the quantity over 
all such non-overlapping cubes within the domain. 

4. Results 

We now present the results of the simulations 
carried out. Each simulation was run for one 
sound crossing time, tc — 3.56 x lO^yrs, of the 
simulation domain. 

Figure [T] shows the density distributions in a 
slice aX X = 0.1 pc (i.e. the mid-plane) and times 
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t — 0.2, 0.5 and 1 tc for the mc-5-512 and the 
mhd-5-512 simulations. Both the non-ideal and 
ideal MHD simulations show some anisotropy with 
respect to the projected initial magnetic field di- 
rection - i.e. the (1, 1) direction - with filaments 
both perpendicular and parallel to this direction. 
This is what would be expected as the material 
will flow preferentially along the (1, 1, 1) direction 
and hence we expect the shock fronts to be nor- 
mal to it. This anisotropy becomes more evident 
as time progresses because the kinetic energy de- 
cays to below the energy of the mean magnetic 
field (which is conserved with time) - hence the 
flow evolves from being dominated by kinetic en- 
ergy to dominated by magnetic energy over the 
lifetime of the simulations. Figure [2] contains a 
plot of the same density distribution for hydro-5- 
256 a.t t ~ 0.2 tc for comparison. It can be seen 
that this is qualitatively different to both the ideal 
and non-ideal simulations shown in figure [1] and 
does not display any signs of anisotropy. 

There are obvious qualitative differences be- 
tween the mc-5-512 and mhd-5-512 simulations at 
all times with the density distribution in the ideal 
case containing much more small scale structure 
than the non-ideal case. This results from the abil- 
ity of the magnetic field to move with respect to 
the fluid. This means that it is not compressed 
as much as in the ideal case, and hence the same 
level of small structure is not formed. Since the 
coupling between the magnetic field and the flow 
field is strong, even though it is imperfect, this 
lack of small-scale structure in the magnetic field 
becomes mirrored in the density field. The im- 
plications of this are discussed in more detail in 
section 14.21 Clearly, in the non-ideal simulation 
energy is not as efficiently transported from large 
scales to smaller scales as in the ideal case. 

We now continue our discussion of the results 
by considering a resolution study to demonstrate 
convergence of our numerical results. We then go 
on to discuss the properties of the density, velocity 
and magnetic field distribution resulting from the 
turbulence, and finally the properties of the energy 
decay. 

4.1. Resolution study 

In order to be reasonably confident that the re- 
sults we present in this work are converged we have 
performed a resolution study for the conditions 



Table 1: Definition of the initial conditions used in 
the simulations in this work. 



Simulation 


Mach number^ 


Resolution 


Comment 


mc-5-64 


5 


64^ 




mc-5-128 


5 


128^ 




mc-5-256 


5 


256^ 




mc-5-512 


5 


512^ 




mc-2.5-512 


2.5 


512^ 




mc-10-512 


10 


512^ 




ambi-5-512 


5 


512^ 


rH = 


hall-5-512 


5 


512^ 


rA = 


mhd-5-512 


5 


512^ 


Ideal MHD 


hd-5-256 


5 


2563 


Hydrodynamic 



"Initial rms Mach number of the flow 




Fig. 1. — Logscale plot of a slice in the mass den- 
sity at times t — 0.2, 0.5 and Itc (top to bottom). 
The left hand column is the non-ideal simulation 
(mc-5-512) and the right hand column is the ideal 
MHD simulation (mhd-5-512). 
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used in the simulations. We have run simulations 
which are identical in all respects except for the 
resolution, which varies from 64^ up to 512'^ (see 
table [1]). We then analyze these in an identical 
fashion and consider the differences between our 
results for the different resolutions. 

Figure [3] shows plots of the kinetic energy nor- 
malized to its initial value as a function of time 
for each of the simulations in our resolution study. 
The 64'^ and 128^ plots do appear to be signifi- 
cantly different from the higher resolution simu- 
lations. However, the 256"^ and 512^ simulations 
are much more similar with a maximum difference 
in the total kinetic energy in the simulation at 
any one time being less than 10%. Table [5] con- 
tains least squares fits of the decay over the time- 
interval [0.2tc, tc] assuming it to be of the form 
for each of the kinetic energy (/3k), the mag- 
netic energy (/?b) and their total (/3Tot)- It can be 
seen that (3k varies by approximately 4% over the 
entire range of the resolution study (simulations 
mc-5-64 through mc-5-512). We feel, therefore, 
that we can be fairly confident of the value of this 
exponent when comparing it with the other simu- 
lations presented here. 

The decay of the energy in the magnetic field is 
more rapid than that in the kinetic energy with a 
difference between /3k and /3b of about 4%. At low 
resolution (mc-5-64) /3b is at its highest which is 
a sign that the low resolution is introducing sufh- 
cient numerical viscosity to induce large amounts 
of numerical reconnection. As the resolution is 
increased /3b reduces to approximately 1.36 and 
stays at about this value even up to the maximum 
resolution of 512'^. We can be reasonably confident 
then that resolution is not affecting our estimate 
of the decay rate of the magnetic energy at our 
maximum resolution. 

The decay of the total energy, being derived 
from the decay of kinetic and magnetic energy, is 
also reasonably well converged with a total change 
over the entire range of the resolution study of 
around 4%. 

Note that, while we can be confident from our 
results that the decay rate is converged, we have 
to be more careful when considering our power 
spectra resul ts presented in sectio n s 14.31 and 14.5.21 
For ex ample, 'Kritsuk et al.' (200f) , 'Sch midt et"al 
(|2009l ) and Lcmastcr & Stone (2009.) present 
power spectra in ideal MHD which suggest that at 




Fig. 2. — Logscale plot of a slice in the mass den- 
sity at time t = 0.2 tc for the hydrodynamic sim- 
ulation (hd-5-256). 




o.ooi I ' ' ' — ' — ' ' ' ' ' ' ' ' — ' — I 

(I.J 1 

Time (cmiiinj rime.) 

Fig. 3. — Logscale plot of the kinetic energy (nor- 
malized to its initial value) for each of the simula- 
tions in the resolution study. 

Table 2: The values of the exponent for the kinetic, 
magnetic and total energy decay for the simula- 
tions presented in this work. These exponents are 
calculated by fitting the data over the time inter- 
val [0.2tc,tc]. 



Simulation /3k /3b f^Tot 



mc-5-64 


1 


34 


1 


40 


1 


35 


mc-5-128 


1 


33 


1 


35 


1 


34 


mc-5-256 


1 


37 


1 


36 


1 


36 


mc-5-512 


1 


40 


1 


37 


1 


39 


ambi-5-512 


1 


40 


1 


35 


1 


38 


hall-5-512 


1 


25 


1 


18 


1 


22 


mhd-5-512 


1 


26 


1 


19 


1 


23 


hydro-5-256 


1 


10 






1 


10 


mc-2.5-512 


1 


21 


1 


29 


1 


23 


mc-10-512 


1 


42 


1 


39 


1 


41 
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resolutions of 1024'^ the turbulent inertial range 
is established over at most a decade in k, while 
for resolutions of 512'^ this falls to around half a 
decade. 

4.2. Energy decay 

We now discuss the behavior of the kinetic and 
magnetic energy in our non-ideal simulations and 
compare them with those for our ideal simulation. 

4-2.1. Kinetic energy decay 

Figure 0] contains plots of the kinetic energy as 
a function of time for the mc-5-512, ambi-5-512, 
hall-5-512 and mhd-5-512 simulations. It is ap- 
parent that the behavior of mhd-5-512 is rather 
similar to hall-5-512 and the behavior of ambi- 
5-512 is similar to that of nic-5-512. This indi- 
cates that, at least for the kinetic energy decay, 
the Hall effect has little impact. However, there 
is a marked difference between those simulations 
incorporating ambipolar diffusion (mc-5-512 and 
ambi-5-512) and those which do not (hall-5-512 
and mhd-5-512). Ambipolar diffusion clearly in- 
creases the decay rate of the turbulence. This can 
also be seen from the data presented in table [2] 
where the exponents of the simulations containing 
ambipolar diffusion are greater by about 10% than 
those without. 

In general we expect that turbulence will decay 
more rapidly in systems which have higher vis- 
cosity than those without. Clearly the ambipolar 
diffusion, although it corresponds to a "viscosity" 
in the magnetic field, acts in a similar way to the 
usual viscous forces in a fluid when considering 
this system. This can be seen in figure [T] where 
the density features are more smeared out in the 
non-ideal case as would be expected if we sim- 
ply introduced a viscous term into the momentum 
equations. 

It is important to note here that the non-ideal 
diffusive terms in the induction equation will lead 
to enhanced reconnection in the magnetic field. 
Since we assume an isothermal equation of state 
there is no path by which the energy released by 
reconnection can find its way to the kinetic en- 
ergy of the system. There may be significant dif- 
ferences between the effects of the non-ideal terms 
on the decay rate of turbulence in the isothermal 
and non-isothermal regimes in molecular clouds. 



4.2.2. Magnetic energy growth and decay 

Figure [5] contains plots of the energy of the 
magnetic perturbations induced by the turbulence 
(normalized to the initial kinetic energy) as a func- 
tion of time for the mc-5-512, ambi-5-512, hall-5- 
512 and mhd-5-512 simulations. A similar trend 
is evident in the results for magnetic energy decay 
as that already noted in section [4. 2. II Again, the 
main result is that ambipolar diffusion increases 
the energy decay rate while the Hall effect does 
little to influence it. This result is borne out by 
the data in table [2] in which the decay exponents 
are greater for the simulations incorporating am- 
bipolar diffusion by about 14% than those which 
do not. 

To gain a little more insight into the interplay 
between the magnetic and kinetic energy in the 
system we now focus on simulations mc-5-512 and 
mhd-5-512 at early times. Figure [S] contains plots 
of the magnetic and kinetic energies of these two 
simulations as a function of time. The magnetic 
energy grows initially as the flow converges and 
compresses the magnetic field in regions through- 
out the computational domain. The kinetic energy 
gradually decays during this time. Once shocks 
form, equipartition between the total magnetic en- 
ergy and kinetic energy is reached and the mag- 
netic field begins to decay. At this time the decay 
of the kinetic energy accelerates due to the dissi- 
pative effect of the shocks which have just formed. 
This dependence on shocks forming is confirmed 
when the growth and decay of the magnetic en- 
ergy is examined for simulations mc-2. 5-512 and 
mc-10-512 and compared with mc-5-512 (see fig- 
ure [7]) . Here we can see that the peak in mag- 
netic energy is reached at later times as the ini- 
tial rms Mach number is decreased. The time at 
which shocks begin to form is determined by the 
initial rms Mach number (giving us a character- 
istic velocity) and half the shortest wavelength in 
the initial velocity field. The approximate shock 
formation times are 0.0125 ic, 0.025 tc and 0.05 ic 
for the mc-10-512, mc-5-512 and mc-2. 5-512 simu- 
lations respectively. These times match up reason- 
ably well, particularly given the temporal granu- 
larity of the simulation data, with the peak of the 
fluctuating part of the magnetic energy in these 
simulations. 

Just after equipartition is reached the decay of 
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Fig. 4. — Logscale plot of the kinetic energy (nor- 
malized to its initial value) for simulations mc-5- 
512, ambi-5-512, hall-5-512 and mlid-5-512. Note 
that the data for the mc-5-512 and ambi-5-512 are 
almost identical. 




Fig. 5. - Logscale plot of the magnetic energy 
(normalized to the initial kinetic energy) for simu- 
lations mc-5-512, ambi-5-512, hall-5-512 and mhd- 
5-512. Note that the data for mc-5-512 and ambi- 
5-512 are almost identical. See text. 



the kinetic energies of mc-5-512 and mhd-5-512 
begin to behave differently. This is unsurprising 
as it is only when the magnetic field perturba- 
tions arc reasonably strong that non-ideal effects 
can play a dynamically significant role. Just af- 
ter the magnetic energy peaks the kinetic energy 
of simulation mhd-5-512 goes through a short pe- 
riod during which it does not decay particularly 
rapidly - presumably a result of the transfer of en- 
ergy from the compressed magnetic field back to 
the kinetic energy as the fluid expands after the 
initial compressions in the simulation. This does 
not happen in simulation mc-5-512 since the mag- 
netic field does not become so compressed in the 
first place due to non-ideal effects. This can be 
seen from the fact that the magnetic energy peaks 
at a lower level in this simulation. In addition, 
energy from this field will not be transferred effi- 
ciently back to the kinetic energy because of the 
imperfect coupling between the magnetic and ve- 
locity fields resulting from the non-ideal terms in 
the induction equation. 

The picture which emerges is the following: ini- 
tially kinetic energy is transferred into magnetic 
energy through compressions. If there are diffu- 
sive terms present in the induction equation then 
this extra magnetic energy will be dissipated and 
there will be less energy available to transfer back 
to the kinetic energy as the compressed regions 
expand. This process will be repeated through- 
out the simulation as shocks form and dissipate. 
Hence we expect a more rapid decay of both mag- 
netic and kinetic energy when diffusive terms are 




Time (crossing times) 

Time (crossing times) 



Fig. 6. — Logscale plot of the magnetic and ki- 
netic energies (normalized to the initial kinetic en- 
ergy) for simulations mc-5-512 and mhd-5-512 for 
a larger range of times. See text. 



Fig. 7. — Plot of the growth and decay of magnetic 
energy (normalized to the initial kinetic energy) 
for simulations mc-10-512, mc-5-512 and mc-2.5- 
512. 
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present. Recalling that the Hall effect does not 
actually diffuse magnetic energy, it is no surprise 
that simulation hall-5-512 is so similar in terms of 
energy decay to mhd-5-512. 

Finally, we note that the peak of the fluctuating 
part of the magnetic energy occurs at a slightly 
later time in the mhd-5-512 simulation than the 
mc-5-512 one. We explain this as follows. As 
noted in Section [4. 3. 11 the ambipolar diffusion ap- 
pears to set a dissipation length scale correspond- 
ing to k — 10. This is significantly larger than 
the dissipation scale in the ideal mhd simulation. 
Hence energy in the magnetic field need not cas- 
cade to such short length scales in order to be 
dissipated in the mc-5-512 and ambi-5-512 simu- 
lations as in the mhd-5-512 and hall-5-512 simula- 
tions. Since the energy cascade to shorter length 
scales takes time (particularly at very early times 
when shocks are only beginning to form), the en- 
ergy in the magnetic field can be removed from the 
former simulations at earlier times than from the 
latter simulations as it need not cascade so far. 

4.3. Power spectra 

We now turn to power spectra of the simula- 
tions. These tell us about the scale of the struc- 
tures formed and are useful in gaining some insight 
into the effects of the non-ideal physics incorpo- 
rated. It is important to stress that since these 
simulations are of decaying turbulence, rather 
than driven turbulence, we do not necessarily 
expect to get the often-quoted power-law depen- 
dence of power on wavenumber. All the power 
spectra presented have been calculated at t — tc 
when we can be fairly confident that the turbu- 
lence is fully developed and the initial conditions 
have been effectively forgotten. In addition, the 
concerns noted in section 14.11 should be borne in 
mind - i.e. the turbulent inertial range is likely 
to be of order half a decade in k for our 512'^ 
simulations. 

We discuss the power spectra of density, mag- 
netic field strength and velocity in turn. 

4-. 3.1. Density power spectra 

Figure [8] contains plots of the spherically in- 
tegrated power spectrum of the density for the 
mc-5-512, hall-5-512, ambi-5-512 and mhd-5-512 
simulations at i = tc. It is clear that there is 



little difference between the mc-5-512 and ambi- 
5-512 simulations - again strongly indicating that 
the Hall effect has little influence on the behavior 
of the simulations. The mhd-5-512 and hall-5-512 
simulations are rather similar, although there is 
less power in small-scale structures in the hall-5- 
512 simulation than the mhd-5-512 one. We leave 
discussion of this until section 14.3.31 

Even at large length scales there is a signif- 
icant difference in the power spectra, indicating 
that dissipation on short length scales due to the 
non-ideal terms in the induction equation indeed 
affects the behavior of the density distribution at 
large scales. The simulations with ambipolar dif- 
fusion have considerably steeper slopes than those 
without (see Tabl e [H in qualita tive agreement 



with the results of lLiet all (|2008h . That the mc- 



5-512 simulation has less power at shorter length 
scales than mhd-5-512 can be seen by a cursory 
examination of figure [1]- the mc-5-512 images ap- 
pear more "blurry" than those from mhd-5-512 - 
so this result is not a surprise. 

Further structure in the power spectrum for 
mc-5-512 is apparent which is absent in the mhd- 
5-512 simulation: at A; « 10 there is a break to 
a steeper slope in the simulations containing am- 
bipolar diffusion. For the mhd-5-512 simulation, 
at about A; « 35 numerical viscosity begins to af- 
fect the power spectrum, as evidenced by the roll 
over of the spectrum from a power-law to a steeper 
slope. The same effect does not appear in the 
mc-5-512 simulation because the physical viscos- 
ity due to the non-ideal terms in the induction 




Fig. 8. — Logscale plot of the normalized spheri- 
cally integrated power spectrum of the density dis- 
tribution for the mc-5-512, hall-5-512, ambi-5-512 
and mhd-5-512 simulations at t = tc. 
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equation dominates the numerical viscosity up to 
much higher values of k. While the power spec- 
trum for the latter simulation does change slope 
at /c ~ 10 it then maintains a strong power law up 
tok^ 100. 

Hence the break at fc « 10 for mc-5-512 
and ambi-5-512 appears to be a physical re- 
sult, rather than a numerical one. This suggests 
that, apparently in cont radiction with a result of 
Oishi fc Mac Low ( 20061 ). ambipolar diffusion can 
set a length scale in turbulence - in these sim- 
ulations that length scale is about 0.02 pc. We 
discuss this further in section 14.3.31 but note that 
we should be cautious about drawing general in- 
ferences given that the resistivities used in these 
simulations are constant in space and time. 

Ultimately, at very short length scales (fc > 100 
or / < 0.002 pc) all four simulations steepen signif- 
icantly. This is expected since at these values of 
fc, which correspond to lengths of less than about 
5 grid zones in these simulations, numerical vis- 
cosity will certainly dominate structure genera- 
tion/dissipation. 



Table 3: The values of the exponent for the power 
spectra of density, velocity and magnetic field 
measured at t = tc- AH fits are over the range 
5 < fc < 20 unless otherwise noted. 



Simulation 


Density 


Velocity 


Magnetic field 


mc-5-512 


2.09^ , 4.06'^ 


1.47 


2.17*^ , 4.96*^ 


ambi-5-512 


2.04^ , 4.03*^ 


1.49 


2.21^ , 4.75*^ 


hall-5-512 


1.41 


1.20 


1.65 


mhd-5-512 


1.45 


1.17 


1.59 


mc-2.5-512 


2.37'^ , 3.32*^ 


2.14 


2.86^ , 3.80^ 


mc-10-512 


2.20'^ , 4.20*^ 


1.44 


2.59^^ , 4.70^ 



"Fitted over 4 < A; < 10 
''Fitted over 10 < fc < 100 



4-3.2. Velocity power spectra 

We now turn to the velocity power spectra. Fig- 
ure [9] contains plots of the spherically integrated 
power spectrum of the velocity for the mc-5-512, 
hall-5-512, ambi-5-512 and mhd-5-512 simulations 
at t = tc. In broad terms of the differences be- 
tween the simulations we obtain the same result 
as for the density power spectrum (cf figure [8|). It 
is, however, obvious that the density and velocity 
power spectra differ quite substantially in qualita- 
tive terms. 

At very low fc the velocity spectra have a rather 
shallow slope with little difference between any of 
the simulations up to fc w 4 - the maximum value 
of fc at which the initial velocity field was non- 
zero. In the range 5 < fc < 20 both mhd-5-512 
and hall-5-512 follow a similar power law, while 
mc-5-512 and ambi-5-512 attain a steeper slope 
in this range. In general, for fc > 5 the simula- 
tions with ambipolar diffusion have significantly 
steeper power spectra, as observed in the density 
power spectra and in qualitati ve agreernent w ith 



( 200ah . 



the driven turbulence results of iLi et al. 

There is a further apparent break in the power 
spectrum at around fc « 100 in the simulations 




Fig. 9. — Logscale plot of the normalized spher- 
ically integrated power spectrum of the velocity 
distribution for the mc-5-512, hall-5-512, ambi-5- 
512 and mhd-5-512 simulations at t = tn. 



10 



with ambipolar diffusion occurring where numeri- 
cal viscosity begins to dominate the ambipolar dif- 
fusion (see section |4.3.1|) . It is interesting to note 
that the series of breaks in the velocity power spec- 
tra in mc-5-512 and ambi-5-512 are not mirrored 
in the density power spectra indicating a level of 
decoupling between these two fields. 



iKritsuk et al.l (|2007f l noted that under certain 
circumstances a power spectrum of p^^^q could 
follow a power law with the classical Kolmogorov 
slope of -5/3 even for compressible (hydrody- 
namic) turbulence. Figure [10] contains power 
spectra of this variable for the mc-5-512 and mhd- 
5-512 simulations. It is apparent that the power 
spectra for this variable for these simulations is 
shallower than the -5/3 law, at least up to the 
value of k where dissipative effects may be im- 
portant. In the range 4 < /c < 10 the exponents 
are -1.289 and -1.14 for the mc-5-512 and mhd-5- 
512 simulations, respectively. The theory under 
which these variables could be expected to have 
the Kolmogorov slope has the assumption that 
the system has reached a statistical steady state - 
since we investigate decaying turbulence here this 
is unlikely to be the case. It might be expected 
that if we do not continually supply energy at low 
k then the slope after a turbulent crossing time 
when the turbulence is well established and has 
also decayed significantly would be too shallow as 
the power spectra would be "too low" at low k. 




Fig. 10. — Logscale plot of the normalized spher- 
ically integrated power spectrum of p^^^q for the 
mc-5-512 and mhd-5-512 simulations at i = tc- 
Also shown is the fc"'^/'^ Kolmogorov power-law. 



4-3.3. Magnetic field power spectra 

Figure [TT] contains plots of the power spectrum 
of the magnetic field at time t = t^. Once again 
these plots indicate that ambipolar diffusion is the 
dominant non-ideal effect in molecular cloud tur- 
bulence. We can see that there are some differ- 
ences between each of the simulations. The mc- 
5-512 and ambi-5-512 simulations are rather sim- 
ilar and are rather similar to the density power 
spectra (figure [8]). However, at high k the ambi- 
5-512 simulation does have noticeably more power 
than mc-5-512. This indicates that the Hall ef- 
fect, present in the mc-5-512 simulation, is having 
some impact in the structure of the magnetic field 
at short length scales. This effect is very small 
in the power spectra of both velocity and density 
and is not apparent in the energy decay rates ei- 
ther. It would appear that when the ambipolar 
diffusion is strong in comparison to the Hall ef- 
fect, as is the case for molecular clouds, the Hall 
effect may change the structure of the magnetic 
field but this change will not propagate into the 
rest of the fluid variables. 

The power spectra for mhd-5-512 and hall-5- 
512 are rather similar to each other, again with 
the exception of high k where the hall-5-512 simu- 
lation has significantly less power. We explain this 
as follows. While the Hall effect does not actually 
diffuse the magnetic field it does re-orient it. This 
process of re-orientation can give rise to a topol- 
ogy favoring some magnetic reconnection - partic- 
ularly on the small scales at which the Hall effect 
operates. We therefore attribute the lower power 
at high k in hall-5-512 to the interplay between 
the Hall effect and numerical viscosity at these 
length scales. Recall that since hall-5-512 contains 
no ambipolar diffusion and negligible parallel re- 
sistivity any reconnection which occurs must be 
almost entirely due to numerical viscosity. The 
same argument holds for the difference between 
mc-5-512 and ambi-5-512 at high fc, but in this 
case the interplay is between the Hall effect and 
ambipolar diffusion. The marked turn up in power 
for k > 200 for the mhd-5-512 simulation is due to 
numerical effects. 

Overall the plots are rather similar to the den- 
sity power spectra (figure[8]): the density and mag- 
netic field distributions appear to be fairly well 
coupled and different in nature to the velocity dis- 
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tribution. 



4.4. Velocity dispersion 

It has been generally accepted that observations 
of line of sight velocity dispersion in molecular 
clouds exhibits a power law with the size of the 
field of view (e.g. lLars"onlll981 ). It is worth noting, 
however, that recent observational re sults call the 



so-ca lled Larson's law into question (jHever et al 
20091) . However, for completeness we feel it is 



worthwhile to examine how our non-ideal simu- 
lations behave in this regard. In this section we 
examine the velocity dispersion as a function of 
characteristic length for each of the simulations. 

Figure [12] contains plots of the velocity dis- 
persions a,t t — tc for each of the 512'^ simula- 
tions. Once again, the results follow the general 
trend of indicating that the Hall effect has al- 
most no impact with the data for mhd-5-512 and 
hall-5-512 being almost indistinguishable. At all 
length scales the velocity dispersion for the sim- 
ulations containing ambipolar diffusion (mc-5-512 
and ambi-5-512) is lower, with the difference being 
larger at short length scales. 

From the power spectra study in section [4?3] we 
already know that the power in short scale varia- 
tions of the velocity is decreased by the presence 
of ambipolar diffusion. It is therefore no surprise 
that comparison of the results for the velocity dis- 
persion shows that the dispersion is decreased by 
the presence of this diffusion. 

There is no obvious power law in these results, 
in agreement with t he simulations of d r iven, ideal 
MHD turbulence of iLemaster fc Stond (|2009f) . As 
noted by these authors, the varying strength and 
direction of the magnetic field give rise to a large 
number of signal speeds within the domain, rather 
than the simple sound speed frequently assumed 
in hydrodynamic turbulence and hence there is no 
reason to expect a power law in the velocity dis- 
persion. 

4.5. The effect of initial rms Mach number 

Finally we turn to the issue of the initial Mach 
number chosen in the simulations. To study this 
we consider simulations mc-2. 5-512, mc-5-512 and 
mc-10-512. These correspond to initial sonic Mach 
numbers of 2.5, 5 and 10, respectively and Alfvenic 
Mach numbers of 0.96, 1.9 and 3.85. 




Fig. 11. — Logscale plot of the normalized spher- 
ically integrated power spectrum of the magnetic 
field strength distribution for the mc-5-512, hall- 
5-512, ambi-5-512 and mhd-5-512 simulations at 

t = U,. 




Fig. 12. — Logscale plot of the velocity disper- 
sion for the mc-5-512, hall-5-512, ambi-5-512 and 
mhd-5-512 simulations. Note that the data for 
simulations mhd-5-512 and hall-5-512 are almost 
identical. 
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4-. 5.1. Energy decay 



Figure [13] contains plots of the decay of kinetic 
energy for each simulation. The energy in the 
plots has been normalized to the starting energy 
for each simulation. It is clear that the higher 
Mach number flows lose their energy more rapidly 
than the lower Mach number ones. This behavior 
is repeated in the magnetic energy and the total 
energy. Table [5] contains the power law indices of 
decay for the kinetic, magnetic and total energy 
in these simulations. The exponents for the ki- 
netic energy decay for the mc-2. 5-512, mc-5-512 
and mc-10-512 simulations are 1.21, 1.40 and 1.42 
respectively, while those for magnetic energy are 
1.29, 1.37 and 1.39 respectively. It is clear that, in- 
deed, all the indices increase with initial rms Mach 
number. We explain this by noting that high Mach 
number flows tend to have strong shocks which will 
dissipate energy more effectively than lower Mach 
number flows. 

4.5.2. Power spectra 

Figures [141 [El ^^^d \W\ contain the power spec- 
tra for the mc-10-512, mc-5-512 and mc-2. 5-512 
simulations for each of the density, velocity and 
magnetic field respectively taken at t — tc. 

The density power spectrum does not appear to 
change much between the mc-10-512 and mc-5-512 
simulations, although it is markedly shallower for 
the mc-2. 5-512 simulation. It is worth noting that 
the mc-2. 5-512 simulation is initially sub-Alfvenic 
and so is qualitatively different in nature to the 
other two. In addition, the rms sonic Mach num- 
ber drops below 1 at f « 0.1 He and hence we have 




Fig. 13. — Logscale plot of the normalized kinetic 
energy for the mc-10-512, mc-5-512 and mc-2. 5- 
512 simulations. 



Fig. 14. — Logscale plot of the power spectrum of 
the density for the mc-10-512, mc-5-512 and mc- 
2.5-512 simulations. 




Fig. 15. — Logscale plot of the power spectrum 
of the velocity for the mc-10-512, mc-5-512 and 
mc-2. 5-512 simulations. 




Fig. 16. — Logscale plot of the power spectrum 
of the magnetic field for the mc-10-512, mc-5-512 
and mc-2. 5-512 simulations. 
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rather well-evolved subsonic, decaying turbulence 
in this case. 

The velocity power spectrum for the mc-10-512 
and inc-5-512 simulations are also remarkably sim- 
ilar, while the mc-2. 5-512 simulation is steeper at 
low k and shallower for 20 < k < 100. Once again, 
we attribute this difference in behavior to the fact 
that the initial rnis velocity of the mc-2. 5-512 sim- 
ulation is sub-Alfvenic. 

The magnetic power spectrum follows a similar 
pattern in terms of differences between the mc-10- 
512, mc-5-512 and mc-2. 5-512 simulations as for 
the other two sets of power spectra. 

4-5.3. Velocity dispersion 

Figure [IT] shows plots of the velocity dispersion 
at t — tc for the mc-10-512, mc-5-512 and mc-2. 5- 
512 simulations. Again, as noted in section 14.41 
no overall power-law is observed. At short length 
scales (less than about 0.01 pc) the slope of the 
relations are all approximately the same. Above 
this scale the slope is somewhat lower for mc-10- 
512 and mc-5-512 than for mc-2. 5-512. This is 
to be expected since large velocity variations will 
be preferentially suppressed by strong shocks dur- 
ing the early evolution of the system, leading to 
each of these simulations having more similar ve- 
locity dispersion at large length scales. Hence we 
expect higher Mach number simulations to retain 
somewhat higher velocity dispersions at all length 
scales than their lower Mach number counterparts, 
but that the fractional differences in these disper- 
sions will be lower at large length scales than at 
shorter ones. 



mf-10-5l2 — I — 
mf-5-512 




Fig. 17. — Logscale plot of the velocity dispersion 
as a function of spatial scale at t = for the mc- 
10-512, mc-5-512 and mc-2. 5-512 simulations. 



5. Conclusions 

We have presented the first study of turbulent 
decay in the presence of non-ideal terms in the in- 
duction equation and also the first simulations to 
incorporate both the Hall effect and ambipolar dif- 
fusion simultaneously in simulations of molecular 
cloud turbulence. This is the first stage in a com- 
prehensive study of non-ideal MHD turbulence. 
The non-ideal effects are therefore included in a 
simplistic way in this work with a view to devel- 
oping a complete, intuitive understanding of their 
impact on turbulence as we continue the study 
and gradually add in more realistic components 
of these effects. The results of this further stage 
of the study are the subject of a forthcoming paper 
currently in preparation. 

We have used a resolution study to determine 
that our simulations are well-resolved for the pur- 
poses of the results which we present. In partic- 
ular, the turbulent energy decay is well resolved 
from a resolution of about 256'^ and the results we 
present here are derived from 512"^ simulations. 

We have found that the Hall effect has little 
influence on kinetic energy decay when present at 
the levels found in molecular clouds. However, am- 
bipolar diffusion increases the rate of energy decay 
on length scales of 0.2 pc or less. Similar results are 
found for the behavior of magnetic energy. Am- 
bipolar diffusion increases the rate of energy decay 
through diffusing away energy while it is stored in 
the magnetic field. Viewed in this way it is un- 
surprising that the Hall effect does not influence 
kinetic energy decay strongly since it does not ac- 
tually diffuse magnetic energy. The kinetic and 
magnetic energy decay is faster in higher Mach 
number flows as would be expected both from this 
argument and from the fact that we are simulat- 
ing isothermal flows which will lose energy more 
efficiently when stronger shocks are present. 

The power spectra of the density for these sim- 
ulations again suggest that the Hall effect has lit- 
tle impact on the flows, while ambipolar diffusion 
cannot be ignored. As might be expected from a 
diffusion term, the power spectrum is softer (i.e. 
steeper) when it is present. At a resolution of 512'^ 
and an assumed length-scale of 0.2 pc we appear 
to have resolved the length at which ambipolar 
diffusion begins to influence the flow. Ambipo- 
lar diffusion does appear to impose a characteris- 
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tic length-scale on the turbulence. However, we 
must be cautious about this interpretation since 
our assumption of spatially constant resistivities 
will have some impact on this result. 

When comparing the velocity power spectra 
with the density power spectra we find that there 
appears to be a decoupling between the two fields 
with breaks in the power-laws not mirrored be- 
tween the two sets of spectra when non-ideal ef- 
fects are included. 

The Hall effect does have some impact on the 
magnetic power spectra at high k. It decreases the 
energy at high k, probably due to re-orientation of 
the magnetic field at small scales in such a way as 
to favor reconnection (either numerical or physi- 
cal) and hence destruction of structure on these 
scales. 

Calculations of the velocity dispersion as a func- 
tion of length scale show that, again, ambipolar 
diffusion is the dominant diffusive term and that it 
has a significant impact. It preferentially reduces 
the velocity dispersion at small scales. We do 
not find a power law dependence between length 
scale and velocity disper sion, in agreement with 
Lemaster fc Stond (|2009t ). 

Comparisons of decaying turbulence with vary- 
ing initial rms Mach number show that higher 
Mach number flows decay more quickly than their 
low Mach number counterparts. There are also 
differences in the power spectra with the Mach 2.5 
flow being significantly different to the Mach 5 and 
Mach 10 flows. The Mach 2.5 flow is slightly sub- 
Alfvenic initially and this may explain the quali- 
tative difference seen. The velocity dispersions are 
lower at short scales for lower Mach number flows, 
but are similar at large scales. This is due to the 
tendency of strong shocks to decay very quickly 
and hence we do not expect even high Mach num- 
ber flows to maintain high velocity dispersion at 
large length scales for long. 

The next step in this work is to study turbulent 
decay in the presence of resistivities determined 
consistently from the magnetic fleld and the den- 
sity of charged species throughout the computa- 
tional domain. This will allow us to understand 
precisely the impact of spatially varying resistiv- 
ities on turbulent decay. We can expect, for ex- 
ample, that the behavior of the power spectra and 
velocity dispersions will be strongly affected. In 



addition, such simulations will then incorporate 
all the non-ideal effects likely to be of importance 
in the study of molecular cloud turbulence. 
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